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Abstract: We measure the energy levels of the excitations of the flux tube between static 
quark and antiquark in three-dimensional SU(2) gauge theory. Combining exponential 
error reduction techniques and a variational method we are able to reduce the errors for 
the excited states significantly and to extract excited states in distinct parity and charge 
conjugation channels. It is conjectured that the infrared behavior (at large qq separation 
R) of the flux tube is governed by an effective string theory. Indeed previous simulations 
show good agreement between lattice data and predictions from Nambu-Goto string theory. 
Recently, new results on the effective string theory obtained corrections to the Nambu-Goto 
predictions and showed that for the open string in three dimensions first corrections should 
appears at order 1/R 4 . They correspond to boundary terms in the worldsheet field theory. 
These corrections are presumably small for the ground state, but significantly larger for 
the excited states and lift the degeneracies of the free theory. Assuming this functional 
form of the correction, we obtain for the coefficient 62 = —0.5(2) (2). 
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1 Introduction 

The formation of a flux tube between static quark and antiquark in the QCD vacuum, 
leads to a linearly rising potential and is one of the possible mechanisms to describe quark 
confinement. Simulations over the last years strongly support this picture (for a review see 
[1]). It is conjectured that at large distances R the dynamics of the flux tube is governed 
by an effective string theory, describing the properties of strong interactions in the infrared 
limit. The construction of models for these QCD strings in terms of open bosonic string 
theories have been attempted for a long time [2]. Today there are two different possible 
approaches for consistent effective string theories. The first idea, due to Polchinski and 
Strominger [3] , was to construct an action consisting of all terms invariant under conformal 
transformations that avoids the conformal anomaly in an arbitrary number of dimensions. 
The first analysis of the closed string spectrum to 0(R~ l ) was extended up to and including 
0(R~ 5 ) in [4-12], showing that the spectrum is equivalent to the spectrum of Nambu-Goto 
string theory (NG) [13, 14] 

E» G (R) = aRJl + ^ (n_±(d_2)) , (1.1) 
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at least to that order. 

The other approach is due to Liischer and Weisz [15, 16], using an action consisting 
of all possible terms constructed directly from the transverse degrees of freedom. The 
corresponding coupling constants are fixed through the duality between open and closed 
strings wrapping around a spatial dimension, which is valid in the presence of at least 
one compactified spatial dimension of length R. Recently their work was extended to 
even higher orders in 1/R and to the open string case [17-19], using the SO(l \ d — 1) 
Lorentz symmetry to constrain additional coupling constants (see also [20]). Leading order 
corrections to the closed string spectrum in d > 2 + 1 possibly appear at order 0(R~ 5 ) 
for the excited states and at 0(R~ 7 ) for the ground state. For the energy spectrum of the 
open string the picture is different due to boundary terms in the effective action. Indeed 
it turns out that first corrections to all states may already appear at 0(R~ 4 ), also for 
d = 2 + 1. In addition, there are possible corrections at 0(R~ 5 ) for d > 2 + 1. It is thus 
natural to look for boundary corrections in the d = 2 + 1 theory, since these are the only 
ones contributing up to 0(R~ 6 ) at least. 

The predictions from these effective theories can be tested by comparison to simulations 
of pure Yang-Mills theories and a lot of studies were done using different theories in 2 + 1 
and 3 + 1 dimensions, see [21] and references therein. In this context it is important to 
note, that liftings of degeneracy in the open string energy levels for d = 3 + 1 in pure 
SU(3) gauge theory where already reported in [22-24]. For recent results on the closed 
string spectrum in 2 + 1 and 3 + 1 dimensional SU(N) gauge theory see [25, 26]. In the 
last years accurate measurements for the width of the flux tube became feasible as well and 
results with improved systematics and high accuracy coincide with the NG predictions for 
this quantity [27, 28]. 

It is an ambitious project to compare the energy spectrum of the flux tube with the 
predictions of the effective model. The stringy behavior is expected to set in at relatively 
large distances and the expected deviations from NG energy levels are small. Thus one 
has to measure a small effect at relatively large energies. In addition one has to control 
systematic effects as e.g. contributions from excited states and effects coming from finite 
extent of the lattice. Indeed most simulations so far haven't been able to reduce the errors 
and to control the systematic effects sufficiently at the same time. This is especially true 
for open strings with Dirichlet boundary conditions, corresponding to the excited states of 
a flux tube between static quark and antiquark. The corresponding energy spectrum can 
be explored using Wilson loops, but large temporal extent is needed to ensure a sufficient 
suppression of contaminations from excited states. 

In [29, 30] a new method was proposed to measure large loops with high accuracy, based 
on the Liischer- Weisz multilevel algorithm [31], and has been applied to the spectrum of the 
flux tube in three-dimensional SU(2) gauge theory. In this study this method is combined 
with a variational method to extract excited state energies in a specific (C, P)-channel 
and reduce excited state contaminations as much as possible. This enables us to extract 
the energy levels accurate enough to compare to the predictions from the effective theory. 
Again we use SU{2) gauge theory and work in three dimensions, but focus only on one 
lattice spacing at (3 = 5.0 with a relatively large lattice spacing in terms of ro, the Sommer 
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scale [32]. Simulating three instead of four dimensions is convenient because simulations 
are much faster and 1/i? 5 corrections to the NG energy levels are absent in the three- 
dimensional theory. 

The paper is organised as follows: In the next section we begin by summarising the 
corrections to the Nambu-Goto energy levels and lay the groundwork for the comparison 
to the measurements. In section 3 we turn to the details of the simulations, describe the 
extraction of the excited states and the control of systematic errors. In section 4 we discuss 
our results and compare to the NG predictions. In section 5 we compare to the boundary 
corrections and finally draw our conclusions in the last section. 

2 Corrections to Nambu-Goto from boundary terms 

In this section we discuss the boundary corrections to the energy spectrum of open strings in 
d = 2 + 1. To keep things simple, we only discuss the main points and adopt the notation 
from [18, 19]. As shown there, the only additional contribution to NG energy levels to 
0(R~ 5 ) in the effective theory in 2 + 1 dimensions is the contribution of a boundary term, 
leading to the Hamiltonian 

3/00 , _\ 

\ 71=1 / 

The states of the Fock space contributing to the three lowest energy levels, denoted by 
E n i, where n is the energy level in the free theory and i is the index for the degenerate 
states at that level, are given by: 
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Here P is the parity quantum number of the state and C what we call the quantum number 
of charge conjugation (called transverse and longitudinal parity respectively in [26]). To 
be precise: C is the exchange between quark and antiquark, combined with a change in 
the direction of the gluonic flow. Using the commutation relations 

[a n , a- n ] = n5 nm (2.2) 
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it is straightforward to compute the corrections to the NG energies of these states due 
to the Hamiltonian (2.1). The corrections are given by: 

^ {d _ 2) 

. it 3 ( (d-2) 



i? 4 \ 60 

£2,1 = -&2 8 + V 



i? 4 V 60 

62,2 = -&2 ^7 32 + V 



R A V 60 

As we see, already the ground state obtains a 0(R~ 4 ) correction from the boundary 
terms, which is nevertheless strongly suppressed by a factor 1/60 and therefore hardly 
visible in our simulation, except for the unlikely case that 62 is large. Nevertheless it might 
be visible in high accuracy simulations for the ground state, as e.g. performed in [33, 34]. 
The corrections to the excited states are enhanced by a factor 240 or more compared to the 
ground state correction and are therefore much easier to detect, even though it is harder to 
extract these states numerically In addition the degeneracy of the states at n = 2 in the 
free theory is lifted and the magnitude of the corrections to the two formerly degenerate 
states differ by roughly a factor of 4. 

Important for the comparison between measurements and the predictions is the radius 
of convergence of the expansion of the string energy levels in 1/R, which can be estimated 
by x c = 1/A, the radius of convergence for the expansion of yl + A x around x = 0. 
In the effective string theory the radius of convergence corresponds to a critical length 
Rc \[o = yl/^c below which the expansion of the square root ceases to be convergent. 
For closed strings this critical length is relatively large, R c \fa > 3.4 for d = 2 + 1, while 
for the open string we have 

Rc \fc\ n=1 = 2.45 , R c v^| n=2 = 3.51 . (2.4) 

As can be seen from [30, 35, 36], extracting the energies of a flux tube of this length is 
possible, especially for the first excited state. Nevertheless the coefficient 62 is not known 
ab initio and we have to extract it from the data. We discuss the extraction of 62 and the 
results in section 5. 



3 Details of the simulations 

Wilson loops are the observables suited to consider the excited states of the flux tube and 
were already used in [22-24] and [29, 30, 35, 36] to extract the spectrum of the excitations. 
In this study, we combine the use of correlation matrices with the error reduction method 
from [29, 30] to extract the excitations spectrum in the (C, P)-channels from large Wilson 
loops and suppress contributions from other excited states as much as possible. 
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Figure 1. Sets of operators to construct the correlation matrices in different (C, P) channels. 

3.1 Extraction of the energy spectrum 
3.1.1 Extraction of eigenvalues 

To extract the energy spectrum we use 8 different operator sets, shown in figure 1. From 

these sets we first construct (C, P)-projectors, given by the linear combinations (see also 
[30]) 

cH — cl i c2 q3 q-1 

°i i ' °i °i °i (o i \ 

sr = s}-s*-s? + s? [6A) 

Q h Cl Q2 i Q3 Q-l 

°i ~ D i °i "r °i °i ■ 

With these projectors we construct a correlation matrix in each of the (C, P)-channels. 

To obtain the energies one has to extract the eigenvalues in the limit T — > oo. Extract- 
ing these eigenvalues is a delicate topic, since loops with a temporal extent much larger 
than 3 ro are out of reach with todays algorithms. In addition the number of operators N 
used to construct correlation matrices is finite, leading to the problem that even though 
one has extracted the eigenvalues carefully there might be a remaining mixing with states 
belonging to other energy levels. Usually one would use a generalised eigenvalue problem 
(GEVP), as discussed e.g. in [37]. The GEVP has the advantage that the excited states 
contributing to each of the eigenvalues start at n = N + 1 as shown in [38] . In our case it is 
hardly possible to establish this method accurately, since the correlation matrix might be 
ill-conditioned even for the smallest values of T. This is especially an issue in the case of 
the correlation matrices in the (— , — ) and (— , +)-channels, where the signal to noise ratio 
of the larger eigenvalues becomes very small. We therefore diagonalise the correlation ma- 
trices at each value of T separately to obtain the eigenvalues X^ P (R, T) with n = 0, . . . , 7, 
using the QR reduction method (see e.g. [39]) and check explicitly the T dependence of 
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the overlaps between eigenvalues and operators. This method is analogous to the method 
used in [26] , where a variational criterion has been applied at each value of T seperately. 
For all states discussed in this study there is no significant change in the overlaps in the 
range of 4 < T < 12 and we do not expect to see any relicts of this in the data. In contrast 
to the case of the GEVP we now have a remaining mixing with each energy level in the 
channel, even though the mixing between the smallest eigenvalues should be suppressed if 
we have chosen the operators wisely. In [25] and [26] this mixing is minimised by using 
a very large number of basis states. Increasing the number of operators above a limit of 
around 50 is hardly possible in combination with our error reduction algorithm, since the 
needed computing ressources grows drastically. 

3.1.2 Groundstate energies in a given channel 

In this study we mainly focus on the ground states of the (C, P)-channels, except for the 

first excited state in the (+, +)-channel, which is discussed in the next section. Each of 

the eigenvalues corresponding to the ground state in a channel now obeys the spectral 
representation 

oo 

\<Z P {R,T)=Y J P? P {R)<r E?P{R)T , (3-2) 

3=0 

where E^ p are the energies in the (C, P)-channel and j3j the overlaps with the eigenstate. 

If the excited states are suppressed sufficiently, we can extract the corresponding en- 
ergies using the asymptotic behavior 

-ln(A CP (i?,T)) = e£ p (R)T - In (ffi p (R)) , (3.3) 

which is exact in the limit T — > oo. In practice Eq P (R) is obtained by fitting the data to 
(3.3) together with the logarithm as an additional fit parameter. These energies are called 
naive from now on, because contaminations from excited states are neglected. 

Unfortunately the simulations have shown that the suppression of the excited states is 
not sufficient with the loop sizes in reach with the present methods of error reduction. This 
is especially crucial for larger values of R as shown in section 4. It is therefore mandatory to 
take the corrections into account. This can be done by fitting the eigenvalues for different 
temporal extents T a and T b to the leading order formula (see also [29, 30, 35, 36]) 



1 In 



Xk P (R,T b ) 
tf p {R,T a ) 



E C p {R) + -^—cf p {R)e~ sCP W T < 

lb- la (34) 



Here T a < T b , a CP {R) = ff p {R)/ffi p {R) and 5 CP (R) is the energy gap to the first 
excited state in the channel. The energies Eq P (R) are called corrected from now on 
and are obtained as fit parameters together with the a's and <5's from fits to all possible 
combinations of T a and T b . The control of these fits is discussed in more detail in section 
3.3.2. 
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3.1.3 Excited states in a given channel 

For an eigenvalue corresponding to an excited state in a channel the situation is different 
due to possible mixings with states of smaller energy. The only excited state in a channel 
used in this study is the first excited state in the (+, +)-channel on which we focus in this 
section. The only state with smaller energy contributing is the corresponding ground state, 
whose energy Eq + (R) can be extracted with the procedure described above. For the naive 
energy Ef + we can use a fit to the form (3.3), since mixings with other energy levels are 
neglected in this quantity. Nevertheless the extraction of the corrected energy E^ + needs 
some additional discussion. 

A fit to the form (3.4) removes the contaminations from states with larger energy 
values. This remains true for excited states in a channel, but one has to include also the 
corrections coming from states with smaller energy. For the state Ef + we can check the 
contribution of the state E^ + by adding the corresponding term 



T b -T c 



(E++(R)-E++(R)) T a A _ e (E++(R)-E++(R)) (T b -T a )\ (3 ^ 



(i 



to the fit (3.4). Performing a fit using (3.4) together with (3.5) shows, that the contribution 
of Eq + (R) should be well below our statistical errors, since 7 is in all cases a number below 
10~ 4 , so that the total contribution of (3.5) is below 5- 10~ 3 , well below the statistical errors 
for the corrected energy Ef + . The other fit parameters agree well with the parameters 
obtained by a fit to the form (3.4) within the statistical errors. As corrected values for 
Ef + (R) we thus use the results from a fit to the form (3.4) and neglect the contribution 
of the ground state. 

3.1.4 Energy differences 

The extraction of the energy differences can be treated independently from the extraction 
of the total energies. In this way the energy differences might serve as an independent 
check of the asymptotic results for the total energies. Using (3.2) one can obtain for the 
naive energy differences the formula 



In 



( XOP(R, T) \ _ 1e cp {r) _ E CP' {R) -\ T _ ln ( H \ f m 



m,0 



which is similar to (3.3) for the total energies. This formula for naive differences is also 
valid for the differences between excited states in a given (C, P)-channel. 

As before we are mainly interested in the differences between the energies of the ground 
states in different channels. For these the leading order formula for the corresponding 
corrected energy difference is the same as in [30], 
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Table 1. Run parameters of the simulations, ro is the Sommer parameter, t s the temporal extent 
of the sublattices in the LW algorithm, N s the number of updates of the sublattice containing 
the spatial operators and N t the number of updates of the sublattices containing only the time 
transporters. 



where a is a suitable combination of the overlaps and 5 corresponds to the energy gap to 
the next excited state in the channel CP' (which is equivalent to the gap in the channel 
CP to leading order in 1/-R). In addition to the energy differences between the ground 
state energies we are also interested in the energy differences between Ef + and Eq + . In 
this case the two terms corresponding to the mixing between Ef + and E^ + are a term of 
the form (3.5) and a term of the form 

( ~f( R ) \ e (E++(R)-E++(R)-S(R))T a 

T b -T a \Po(R)J (38) 

X (l- e (Et + (R)-E++{R)-6(R)) [T b -T a )^ < jq-4 

The total contribution of these terms is below 5 • 10~ 3 and thus negligible compared to the 
statistical errors. We therefore use a fit to the form (3.7) for that difference, too. 

3.2 Algorithm and simulation parameters 

The need for suppression of the contributions from excited states and the onset of string- 
like behavior in the large R regime demands large Wilson loops, and an efficient algorithm 
with sufficient error reduction for large loops is needed to extract the spectrum with high 
accuracy. The algorithm used in this study is a variation of the Liischer Weisz algorithm 
[31] and discussed in detail in [30]. 

Our simulations were done using the SU(2) Wilson plaquette action in 2+1 dimensions 
at B = 5.0. The Sommer parameter is known e.g. from [29, 35, 36] and listed together with 
the simulation parameters in table 1. In fact the lattice spacing in terms of the Sommer 
parameter is not very small, but a comparison from the results of the four different B values 
from [29] and [30] shows that there is not much movement in the energies expressed in terms 
of ro over the whole range of 5.0 < B < 12.5. We therefore expect that our results are 
relevant for the continuum as well. In contrast to [29, 30, 35, 36], we work with 5 different 
temporal extents instead of 4, which leads to strong improvements in the extraction of the 
asymptotic behavior and the control of contaminations from excited states via (3.4) and 
(3.7). 
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Statistical errors are estimated using the usual binned jackknife method with 50 bins 
for all measurements. We explicitly checked that none of the error estimates varies more 
than a few percent with bin size. 

3.3 Sources of systematic effects 

3.3.1 Finite volume effects 

The first class of sources of systematic errors in addition to the ones discussed in section 
3.1 is due to finite extent L of the lattice. Here periodic boundary conditions are applied 
to the lattice, making it possible for the loops to interact with themselves by around-the- 
world glueball exchanges. In the string picture these exchanges corresponds to handles on 
the world sheet, wrapping around a compactified dimension. The contribution of such a 
handle to the signal of the Wilson loop takes the form 

a(L') exp (-m G L 1 ) , (3.9) 

where L' is the length of the handle, mg the mass of the lightest glueball (measured in 
[40] for SU(2) in d = 2 + 1 and at /3 = 5.0) and a(L') is the overlap with the loop. In 
principle these handles can wrap around any of the directions, but it is immediately clear 
from (3.9) that the main contribution comes from the direction that enables the shortest 
handle. Clearly this is the direction parallel to the spatial direction of the Wilson loop 
where L' = L — R. The main problem thus arises for Wilson loops with large values of R. 
For the case of our study, 2 + 1 dimensional SU(2) with /3 = 5.0, a test for finite volume 
effects with volumes of 24 3 and 48 3 , excited states up to n = 3 and R < 12 already exists 
[30]. Since no finite volume effects were visible in this test, we can conclude that a(L') is 
at most a factor of 0(1) and that we do not expect to suffer from any finite volume effects 
up to energy levels of n = 3 at least. 

3.3.2 Control of the fits 

In the analysis a lot of fitting is done in order to extract the energies and energy differences 
in the limit T — > oo and to take the contaminations from excited states into account. For 
the control of systematic effects it is thus crucial to control these fits in several ways, since 
these are the biggest source for systematic errors in the final data. 

Even though the fits (3.3) and (3.6) include only a linear fit function it might be that 
already these fits give biased results, since the more important points for the asymptotic 
behavior, the points from loops with large T, come with bigger statistical errors. This 
might lead to the problem that the fit is completely determined by the values at small T. 
As a check we plotted the data along with the resulting lines and checked the deviations 
from the line of the points at large T. Where the deviations were not at the percent level 
and below the statistical errors, we discarded the fit. 

The second type of fits, eq. (3.4) and (3.7), are much harder to control, since these 
fits include nonlinear functions of two arguments. In addition the % 2 /d.o.f. of these fits is 
usually very small and does not provide much information on the goodness of the fits. We 
thus apply two checks that were already discussed in [30]. We expect a to be smaller than 
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Figure 2. Cross check of the fits for the corrected energies (3.4) and energy differences (3.7): Left: 
Fits for Eq , R = 5,8 and 11. These fits are expected to give reliable results. Right: Fits for 
E^ , R = 5,8 and 11. These fits are regarded to be unreliable, since the points are clearly not 
following the A = A-linc. We have not plotted the error bars in the right fit, because they are large 
and confuse the picture. 



the ratio of the degeneracies between the next excited state and the state considered and 
5 to be of the order of the energy gap between the two. Wherever this criterion was not 
fulfilled we did not use the fits. The second check is a comparison between the expected 
corrections 
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(3.10) 



(3.11) 



for each R and all combinations of T a and T b . We expect to obtain A = A for each R. The 
plot of A against A for three examples where we expect the fits to work are shown on the 
left of figure 2, while on the right we show three examples where the fits are regarded to 
be unreliable. Similar checks were performed for the energy differences, too [30]. 



4 Results of the simulations 

We list all results, naive and corrected in tables 2 and 3, to enable other groups to use 
the data for comparisons. In this section we compare our data to the full NG energy 
spectrum (1.1) and its truncations in 1/R 2 to leading order (LO), next-to leading order 
(NLO) and next-to- next-to leading order (NNLO). The comparison to corrections from 
boundary terms is postponed to the next section. 

In order to compare to the data we have to extract the string tension a and fix an 
unphysical constant Vq. Since we expect corrections to appear with respect to the full NG 
prediction (as concluded e.g. in [26] and [19]), we fit the ground state data to the form 
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Figure 3. Results for the ground state rescaled such that Eq' =0. The lines are the LO and 
NG predictions as defined by (1.1) in dimensionless quantities, eq. (4.3). 





£ ++ = E 


E+- = E 1 


E+ + = E 2A 


E — E 2 ,2 


R 


naive 


corr 


naive 


corr 


naive 


corr 


naive 


corr 


4 
5 
6 
7 
8 
9 

10 
11 
12 


0.57281(3) 0.5716(1) 
0.67623(4) 0.6756(2) 
0.77831(5) 0.7775(2) 
0.87969(7) 0.8781(3) 
0.98000(8) 0.9779(4) 
1.08047(9) 1.0772(5) 
1.18005(11) 1.1761(6) 
1.28020(13) 1.2749(8) 
1.37937(15) 1.3734(9) 


1.1278(2) 1.118(1) 
1.1643(2) 1.152(1) 
1.2115(2) 1.198(1) 
1.2735(2) 1.255(1) 
1.3387(2) 1.318(2) 
1.4159(2) 1.387(2) 
1.4899(2) 1.460(2) 
1.5766(2) 1.537(2) 
1.6559(3) 1.616(3) 


1.452(2) 1.35(6) 
1.491(1) 1.41(2) 
1.516(1) 1.45(2) 
1.581(1) 1.51(1) 
1.648(1) 1.56(1) 
1.731(2) 1.62(1) 
1.846(2) 1.69(2) 
1.940(3) 1.76(2) 
2.046(4) 1.84(4) 


1.5056(9) 

1.5110(8) 1.477(12) 
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1.6685(7) 1.648(5) 
1.7296(8) 1.709(4) 
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1.8644(9) 1.843(3) 



Table 2. Results for the total energies in lattice units. In addition to the corrected data, which 
are our final results, we also list the naive results. They give upper bounds for future simulations 
and illustrate the effect of contaminations from excited states. Note that also the naive energies 
are already asymptotic results for T — > oo, neglecting the contaminations. 




Since we expect string like behavior to set in at ~ 2 ro [30, 33, 34], we omit the first three 
points in the fit and obtain 

a = 0.0975(2) and V = 0.2148(7) . (4.2) 

This result for a is consistent with the value from [33, 34] within the error bars. The 
error bars are much larger in our case, which is no surprise since our aim is not to obtain 
the ground state energy and the string tension with high accuracy and our setup is not 
tuned to reduce this error in particular. From now on we are going to use the dimensionless 
quantities 

R* = ^R and E* = (E - V ) fyfa (4.3) 

for plotting and comparisons to the predictions. 

We plot the naive and corrected results for the ground state in figure 3. In the plots 
we have rescaled the energies such that En ,LO = n, i.e. 

E* njrsc (R*) = (E* n (R*) ~ R*) ^ + ^ , (4.4) 
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E o - E o + 
= AE W 


E++ 
= AE 20 


E o ~ E o + 


Et + E+- 
= AE 2 i 


R 


naive 


corr 


naive 


corr 


naive 


corr 


naive 


corr 


4 
5 
6 
7 
8 
9 

10 
11 
12 


0.5558(2) 0.547(1) 
0.4879(2) 0.476(2) 
0.4330(2) 0.421(2) 
0.3936(2) 0.377(2) 
0.3584(2) 0.340(2) 
0.3351(2) 0.310(2) 
0.3094(2) 0.283(2) 
0.2956(2) 0.261(3) 
0.2757(2) 0.242(3) 


0.814(1) 0.734(18) 
0.736(1) 0.672(14) 
0.699(2) 0.628(11) 
0.666(1) 0.585(9) 
0.648(2) 0.539(11) 
0.513(16) 
0.483(19) 
0.467(35) 


0.9334(9) 

0.8341(8) 0.800(12) 
0.7506(8) 0.733(7) 
0.6867(7) 0.671(6) 
0.6325(7) 0.617(5) 
0.5684(8) 0.571(5) 
0.5477(8) 0.533(4) 
0.5122(9) 0.499(3) 
0.4828(9) 0.469(3) 


0.319(1) 

0.296(1) 0.25(2) 
0.296(2) 0.25(1) 
0.299(1) 0.24(1) 
0.302(2) 0.23(1) 
0.23(2) 
0.22(2) 
0.22(4) 



Table 3. Results for the energy differences in lattice units. As in table 2 we list corrected and 
naive results. For some differences we were able to obtain corrected but no naive results. This is 
due to the fact that the contaminations from excited states are large and the data points for the fit 
(3.6) do not lie on a straight line, whereas the fit (3.7) still works well. 



to provide visibility of small effects. We see a clear splitting between the naive and the 
corrected results, which increases when we go to larger values of R*. The effect becomes 
even more severe when we go to the excited states, as shown for energy level E± in figure 4 
(top left). We can also see the approximate linear asymptotic behavior of the rescaled naive 
energies when going to large R values, which was seen in [22-24], too. The asymptotic 
behavior of the corrected energies is different and we conclude that the behavior of the 
naive energies is due to contributions from excited states and does not reflect the physical 
asymptotic behavior of the energies for large R. When we look at the predictions for 
the energy gaps from (1.1) and its truncations it is easy to see why the excited state 
contamination is enhanced. The energy gaps to the excited states decrease roughly like 
1/R, thus the damping with T is reduced in eq. (3.2) when going to larger R and more 
excited states contribute to the contamination. In the following we therefore only use 
corrected results for the discussions. 

The results for E^ + and Eq , belonging to the second excited state, are shown in 
figure 4 (top right). It is remarkable that with the new method and 5 different temporal 
extents we were able to reduce the errors for the corrected values of the ground state in 
the (+, — )-channel by a factor of more than two and to obtain corrected results for the 
groundstate in the (— , — )-channel with an accuracy better than that of the results for the 
ground state of the (+, — )-channel in [30]. The error bars of the first excited state in 
the (+, +)-channel are almost an order of magnitude bigger than the error bars from the 
ground state of the (— , — )-channel, but nevertheless we see a 2 — 4 sigma splitting up to 
R* ~ 3. At that point the error of E^ becomes to large to distinguish between the two. 

The energy differences are listed in table 3. They have the advantage that the compar- 
ison to the predictions does not involve the unknown constant Vo from eq. (4.1) and that 
they are more sensitive to subleading properties of the flux tube. In figure 4 we show the 
difference Eq — Eq + (left bottom), corresponding to the free difference AE±o, and the 
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Figure 4. Top: Results for the first (left) and second excited state (right), rescaled such that 
E^' LO = n. Bottom: Results for the energy differences corresponding to Ai?i (left) and AE 2 a 
(right) in the free theory. The data is rescaled such that AE*^° = n — m. 



differences E± — Eq and Eq — Eq (right bottom), corresponding to the free difference 
A£^2o- We see a similar error reduction as for the absolute energy values compared to [30]. 

The qualitative agreement between the data and the NG predictions is remarkable 
down to very small values of R*, where we do not expect the flux tube to have a string 
like shape and where the expansion of the string action in 1/R* is no longer expected to 
converge. These findings are in full agreement with other simulations, e.g. [29, 30, 35, 36]. 
In addition they are consistent with the results from [15, 22-24] for smaller values of R* and 
with the results obtained for closed strings [25, 26]. Even though the data follows the NG 
predictions qualitatively, we see a significant deviation for the absolute energies Eq and 
Eq , as well as for the corresponding differences to the ground state, even for our largest 
values of R* . We are going to compare these deviations to the boundary corrections in the 
next section. The state is fully consistent with the NG prediction inside the error 

bars for all values of R*, which is also true for the corresponding difference Ef + — Eq + . 



5 Discussion of 1/i? 4 corrections 

To compare the boundary corrections to the NG predictions with the data, we have to 
extract the unknown coefficient 62- There is a hint from theory that 62 is negative and 
nonvanishing, coming from a computation of 62 for the case of a string ending on two d- 
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a' 


Vq 


b 2 


7o 


7i 


Vo 


m 


X 2 /d.o.f. 


1 


0.0974(1) 


0.2151(4) 


-0.30(4) 










0.79 


2 


0.0975(2) 


0.2145(7) 


-0.61(19) 


2(3) 


-2(1)-10 3 






0.04 


3 


0.0975(2) 


0.2145(6) 


-0.52(14) 






14(11) 


-11(6)-10 3 


0.05 


4 


0.0975(2) 


0.2146(7) 




9(3) 


12(3)-10 3 




-7(2)-10 4 


0.22 



Table 4. Results for the combined fit of ground state and first excited state data to the form (5.1). 
The different fits arc described in the text. 



branes in confining gauge theories with a weakly curved holographic dual [18]. Corrections 
appears at all energy levels and we demand that b% is consistent with all data above the 
critical distance as defined in section 2. We therefore use the ground state as well as the 
first excited state to extract 62 from a simultaneous fit. We also have to account for the 
possibility that the values for a and Vq are changed due to the ground state correction. 
We thus use them again as free parameters denoted as a' and Vq to distinguish them from 
the values obtained in the last section. 

It is known that, when fitting to a polynomial, the coefficient of the highest power of the 
polynomial in the fit receives a summed contribution from possible higher order terms. This 
spoils the reliability of the highest coefficient of the fit polynomial. 62 suffers from the same 
problem, since the fit function is a nonlinear function in a' combined with a polynomial 
correction. To improve the reliability of the obtained value for 62 we include a correction 
term of higher order to each energy level. There are in general two possibilities at which 
order the next correction term might appear. At 0(R~ e ) if there are additional boundary 
terms in the effective action or at 0(R~ 7 ) if the next correction term is a regular (non- 
boundary) correction. To account for these possibilities we parametrise our fit functions 
for a simultaneous fit to the ground state and the first excited state as: 

(5.1) 

. !~ 23 vr ,241 tt 3 1 1 1 

To check the effect of the higher order terms we perform the following fits: 

1. Set 7i = and rji = 0, use 62 as a fit parameter. 

2. Set r\i = 0, use 62 and ji as fit parameters. 

3. Set 7i = 0, use 62 and rji as fit parameters. 

4. Set 62 = and 770 = 0, use ji and rji as fit parameters. 

In all cases a' and Vq are free parameters as well. The last fit accounts for the possibility 
that the coefficient 62 vanishes identically, which is not yet ruled out completely by theory 
and thus remains a possibility. In the fits we include all points except the one with smallest 
R for the ground state, and all points with R > 7 for the first excited state, since we expect 
the data to be consistent with the effective string theory even for smaller values of R when 
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Figure 5. Comparison between fit 2 from table 4 and the data for the groundstate (left) and the 
first excited state (right). The data is rescaled as in the previous figures. The line marked with 
R* is the line below which the expansion of the corresponding energy level in l/R* ceases to be 
convergent. 



the corrections are included. We exclude the data for the second excited state for the 
following reasons: Only the largest values of R* are above the critical length below which 
the expansion of the energy ceases to be convergent for the second excited state (see eq. 
(2.4)) and one has to include four more fit parameters 72,4 and 772,1 to the fits. Therefore the 
inclusion of the second excited state does not improve the information and the accuracy 
for 62. We do not account for the possibility of correction terms below 0(R~ 4 ), since we 
compare to the effective theory where these terms are explicitly ruled out. 

We list the results of the fits in table 4. In all cases we see good agreement between 
a' and Vq' and a and Vq obtained in the previous section. This shows that the extraction 
of a from the ground state is not sensitive to l/R 4 corrections within the accuracy for the 
ground state data of our simulations. The result for 62 varies between fits 1 — 3, but is 
small and negative in all cases. We expect the results for 62 from the fits 2 and 3, including 
also the higher order terms, to be unbiased. In that case 70 or 770 is practically consistent 
with zero, which is not surprising, since we expect the data for the ground state not to be 
accurate enough to make higher order corrections to a l/R 4 correction visible. 71 and 7/1 
are nonzero but still small compared to the NG coefficient at 0{R" 7 ). Fit 4 seems to work 
as well but a relatively large nonvanishing boundary term at 0(R~ e ) is needed for the fit 
to work accurately. Nevertheless the corresponding x 2 /d.o.f. is a magnitude bigger than 
for the fits 2 and 3 (even though it is still below one). Assuming that 62 does not vanish 
identically we find the data to be consistent with 

62 = -0.5(2)(2) , (5.2) 

where the first error is statistical and the second reflects the systematical uncertainty, 
estimated from the variation of 62 between fits 1-3. 

To compare to the data, we use fit 2 from table 4 and show the corresponding plots for 
the ground state and the first excited state in figure 5. We use a and Vq as before (instead 
of a' and Vq) to rescale the data points and show also the NG lines from the previous 
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Figure 6. Comparison between the data for the second excited state and the corrected energy 
levels _E 2 ,i = E^ + £2,1 and E 2l 2 = E^ G + £2,2 with b 2 obtained from fit 2. The data is rescaled 
as in the previous figures. 

section for comparison. It is remarkable that we can accurately describe the corrections to 
two different energy states with a predicted splitting and a single coefficient. In addition 
we also show the resulting curve from fit 4 for the first excited state. We see that it deviates 
a little bit more from the data than the curve obtained from fit 2, but is well inside the 
error bars and describes the data as well. 

Having fixed 62 through the fit (in this case fit 2), we now get a prediction for the two 
curves corresponding to the second excited state, where the degeneracy is lifted. We show 
the corresponding curves together with the data in figure 6. Additional terms of 0(R~ e ) 
or higher are excluded here, since we do not know anything about their coefficients and we 
expect the effects to be of minor importance at the present level of accuracy of the data. 
We see that the curve £2,1 stays below the £2,2 curve which is consistent with the behavior 
of the data. Nevertheless the data is closer to the NG curve when R* < 3, which is well 
below R* for the second excited state. Above R* the data for Eq approaches the £2,2 
curve, but it is not clear whether the behavior with increasing R* remains consistent with 
the resulting curve. The data for Ef + is consistent with the £2,1 curve already below R*, 
but the error bars above R* are to large to exclude the approach to any of the three curves. 

The energy differences provide another check for the fits of this section. We show the 
curves obtained with the parameters from the fit above and the corresponding data in fig. 
7. Above R* we see perfect agreement between the AE'io curve including the correction 
and the data. For AE^o the picture is similar as for the total energies at n = 2, the data 
is roughly consistent with the predictions. 

We can employ yet another independent check and use the splitting between first 
excited state and ground state directly to obtain 62- In that case we can keep o and Vb 
fixed, which are almost completely determined by the groundstate. When we include an 
additional correction term at 0(R~ 6 ) we obtain 62 = —0.56(23) with x 2 /d.o.f. = 0.03. This 
result is consistent with the results from the fits to the total energies. All in all, together 
with the hint from the effective theory that 62 should be nonvanishing, we consider the 
possibility of a vanishing 1 /R 4 correction to be unlikely with our present data. The data 
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Figure 7. Comparison between the data for the energy differences AEio (left) and AE20 (right) 
and the predictions with 62 obtained from fit 2. The data is rescaled as in the previous figures. We 
have used the abbreviations Aeio = ei — £0j Ae2o,i = £ 2,i — £o and Ae20.2 = £2,2 — eo- 

for ground state and first excited state is well described by the predicted splitting of the 
effective theory with a single coefficient 62. Also the behavior of the second excited state, 
which was not included in the fits, suits to the predictions. 

6 Conclusions 

In this article we have looked at the energies of the excited states of the flux tube between 
static quark and antiquark in three-dimensional SU(2) gauge theory, using a combination 
of the algorithm from [30] and a variational method. We studied qq separations between 
1 and 3 ro and energy levels up to n = 2. We have discussed the main systematic effects 
entering the extraction of the excited states in detail, explained under which conditions we 
consider the results to be unaffected by any of these effects and included only those results 
in the analysis where systematics are under control. The only remaining systematic effect 
which is not investigated is the lattice spacing dependence, but results from [29, 30, 35, 36] 
suggests that the dependence on the lattice spacing is very mild. Using the correlation 
matrices in the (C, P) channels and 5 different temporal extents, we were able to further 
reduce the errors compared to [30] and to obtain a clear signal for the second excited state. 
In addition we were able to obtain results for the first excited state in the (+, +) channel, 
which is degenerate to the ground state in the (+, — ) channel in the free theory. The 
effective string theory predicts a splitting of these formerly degenerate states and indeed 
we observe such a splitting in our data. In all cases we see the same qualitative agreement 
with the NG predictions as observed for previous measurements in the three-dimensional 
case, e.g. [29, 30, 35, 36]. 

Nevertheless, with enhanced accuracy and improved control over systematic effects we 
see a significant deviation from the NG curves and compare it with the predictions from 
[18, 19]. Our data is in good agreement with the predicted 1/R boundary correction to 
the effective string theory, even though we cannot exclude the possibility that boundary 
corrections do not appear until 0(R~ 6 ). If we assume that 62 does not vanish identically 
we find evidence that 62 is small and negative, 62 = —0.5(2) (2), where the first error is 
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statistical and the second systematic. The resulting curves are in good agreement with 
the data and it is remarkable that we can describe the data accurately with the predicted 
splitting. The negative sign of the result is consistent with the sign of the result for confining 
gauge theories with a weakly curved holographic dual, applied to the case of a string ending 
on two d-branes [18]. This is the first time that significant deviations from NG energy levels 
could be observed for open strings in three dimensions. In the case of closed strings and in 
the finite temperature behavior of the string tension analogous deviations were observed 
in [26] and [41]. 

To confirm the findings of this paper it is desirable to check the consistency with 
earlier results, as e.g. the results for the ground state from [33, 34], and to perform 
simulations at additional lattice spacings. One would also like to have a precise prediction 
for the coefficient 62 to compare with, which is hard to get if 62 is not constrained by any 
symmetry. In addition it would be interesting to see whether 62 is universal or whether it 
varies between different theories. To clarify the picture one also has to push the simulations 
to bigger values of R, in order to reach the region of a convergent expansion for the higher 
energy levels, and to further increase the precision for the excited states. 
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